Theoretical insights into the HO●-induced oxidation of chlorpyrifos pesticide: Mechanism, kinetics, ecotoxicity, and cholinesterase inhibition of degradants.

The oxidation of the common pesticide chlorpyrifos (CPF) initiated by HO● radical and the risks of its degradation products were studied in the gaseous and aqueous phases via computational approaches. Oxidation mechanisms were investigated, including H-, Cl-, CH3- abstraction, HO●-addition, and single electron transfer. In both phases, HO●-addition at the C of the pyridyl ring is the most energetically favorable and spontaneous reaction, followed by H-abstraction reactions at methylene groups (i.e., at H19/H21 in the gas phase and H22/H28 in water). In contrast, other abstractions and electron transfer reactions are unfavorable. However, regarding the kinetics, the significant contribution to the oxidation of CPF is made from H-abstraction channels, mostly at the hydrogens of the methylene groups. CPF can be decomposed in a short time (5-8 h) in the gas phase, and it is more persistent in natural water with a lifetime between 24 days and 66 years, depending on the temperature and HO● concentration. Subsequent oxidation of the essential radical products with other oxidizing reagents, i.e., HO●, NO2●, NO●, and 3O2, gave primary neutral products P1-P15. Acute and chronic toxicity calculations estimate very toxic levels for CPF and two degradation products, P7w and P12w, in aquatic systems. The neurotoxicity of these products was investigated by docking and molecular dynamics. P7w and P12w show the most significant binding scores with acetylcholinesterases, while P8w and P13w are with butyrylcholinesterase enzyme. Finally, molecular dynamics illustrate stable interactions between CPF degradants and cholinesterase enzyme over a 100 ns time frame and determine P7w as the riskiest degradant to the neural developmental system.

Cyclo[n]carbons and catenanes from different perspectives: disentangling the molecular thread.

All-carbon atomic rings, cyclo[n]carbons, have recently attracted vivid attention of experimentalists and theoreticians. Among them, cyclo[18]carbon is the most studied system. In this paper, we summarize and review various properties of cyclo[n]carbons, emphasising the aspects of their aromaticity/antiaromaticity. In the first part, the trends in bonding patterns and selected aromaticity indices with the increasing size of the rings are discussed. In the second part we explore the properties of catenane models based on interlocked cyclo[18]carbon rings from different perspectives and investigate their behaviour under the action of external force using computational experiments.

New insight into environmental oxidation of phosmet insecticide initiated by HO˙ radicals in gas and water - a theoretical study.

Phosmet is an organophosphorus insecticide widely used in agriculture to control a range of insects; recently, it was banned by the European Union in 2022 due to its harmful effects. However, its environmental degradation and fate have not yet been evident. Thus, phosmet oxidation by HO˙ radicals was theoretically studied in this work using the DFT approach at the M06-2X/6-311++G(3df,3pd)//M06-2X/6-31+G(d,p) level of theory. Three different mechanisms were considered, including formal hydrogen transfer (FHT), radical adduct formation (RAF), and single electron transfer (SET). The mechanisms, kinetics, and lifetime were studied in the gas and aqueous phases, in addition to its ecotoxicity evaluation. The results show that FHT reactions were dominant in the gas phase, while RAF was more favourable in the aqueous phase at 298 K, while SET was negligible. The branching ratio indicated that H-abstractions at the methyl and the methylene groups were the most predominant, while the most favourable HO˙-addition was observed at the phosphorus atom of the dithiophosphate group. The overall rate constant values varied from 1.2 × 109 (at 283 K) to 1.40 × 109 M-1 s-1 (at 323 K) in the aqueous phase and from 6.29 × 1010 (at 253 K) to 1.32 × 1010 M-1 s-1 (at 323 K) in the gas phase. The atmospheric lifetime of phosmet is about 6 hours at 287 K, while it can persist from a few seconds to several years depending on the temperature and [HO˙] concentration in the aqueous environment. The QSAR-based ecotoxicity evaluation indicates that phosmet and its degradation products are all dangerous to aquatic organisms, although the products are less toxic than phosmet. However, they are generally developmental toxicants and mutagenicity-negative compounds.

Hydroxyl radical-initiated decomposition of metazachlor herbicide in the gaseous and aqueous phases: Mechanism, kinetics, and toxicity evaluation.

The oxidation of widely-used herbicide metazachlor (MTZ) by hydroxyl radical (HO•) in the gas and the aqueous phases was investigated in terms of mechanistic and kinetic behaviors using the M06-2X/6-311++G (3df, 3pd)//M06-2X/6-31 + G (d,p) level of theory over the temperature range 250-400 K. The formal hydrogen transfer, HO•-addition, and single electron transfer mechanisms were considered. The overall rate constants in the gas phase range from 8.40 × 1010 to 8.31 × 109 M-1 s-1 at the temperature from 250 to 400 K, respectively, while the ones in the aqueous phase are close to diffusion-controlled rates, with diffusion-corrected rate constants being 1.31 × 109 to 1.27 × 109 M-1 s-1. The formal hydrogen transfer mechanism is the most dominant in the gas phase, whereas the HO•-addition is the most favorable in the aqueous phase. The H-abstraction at two methyl groups and the HO•-addition to C11 and C12 atoms (pyrazole ring), C16 and C18 atoms (benzyl ring) are significant. The short lifetime in the environment, equal to only 4.16 h, requires more attention to this herbicide compound, whereas its lifetime in the aqueous condition varies sharply from half second to several thousand days depending on the HO• concentration. The ecotoxicity estimation of MTZ and its principal transformation products to aquatic organisms suggests that they are harmful or toxic substances. Moreover, the MTZ is a developmental toxicant and mutagenicity-positive, while its decomposed products are developmental toxicants with no mutagenic toxicity. Their bioaccumulation in aquatic organisms is negligible.

Together, Not Separately, OH and O3 Oxidize Hg(0) to Hg(II) in the Atmosphere.

Mercury, a highly toxic metal, is emitted to the atmosphere mostly as gaseous Hg(0). Atmospheric Hg(0) enters ecosystems largely through uptake by vegetation, while Hg(II) largely enters ecosystems in oceans and in rainfall. Consequently, the redox chemistry of atmospheric mercury strongly influences its fate and its global biogeochemical cycling. Here we report on the oxidation and reduction of Hg(I) (BrHg and HOHg radicals) in reactions with ozone and how the electronic structure of these Hg(I) species affects the kinetics of these reactions. The oxidation reactions lead to YHgO· + O2 (Y = Br and OH), while the reduction reactions produce Hg(0), OY, and O2. According to our calculations with CCSD(T), NEVPT2, and CAM-B3LYP-D3BJ, the kinetics of both oxidation reactions are very similar and much faster than their reduction counterparts. Past modeling of field data has supported the idea that OH and/or O3 (rather than Br) dominates Hg(II) production in the continental boundary layer. Almost all models invoking OH- and ozone-initiated oxidation of Hg(0) assume that these reactions produce Hg(II) in one step, despite the lack of plausible mechanisms. The two-step mechanism of formation of HOHg followed by its reaction with ozone helps reconcile modeling results with mechanistic insights.

Improved Mechanistic Model of the Atmospheric Redox Chemistry of Mercury.

We present a new chemical mechanism for Hg0/HgI/HgII atmospheric cycling, including recent laboratory and computational data, and implement it in the GEOS-Chem global atmospheric chemistry model for comparison to observations. Our mechanism includes the oxidation of Hg0 by Br and OH, subsequent oxidation of HgI by ozone and radicals, respeciation of HgII in aerosols and cloud droplets, and speciated HgII photolysis in the gas and aqueous phases. The tropospheric Hg lifetime against deposition in the model is 5.5 months, consistent with observational constraints. The model reproduces the observed global surface Hg0 concentrations and HgII wet deposition fluxes. Br and OH make comparable contributions to global net oxidation of Hg0 to HgII. Ozone is the principal HgI oxidant, enabling the efficient oxidation of Hg0 to HgII by OH. BrHgIIOH and HgII(OH)2, the initial HgII products of Hg0 oxidation, respeciate in aerosols and clouds to organic and inorganic complexes, and volatilize to photostable forms. Reduction of HgII to Hg0 takes place largely through photolysis of aqueous HgII-organic complexes. 71% of model HgII deposition is to the oceans. Major uncertainties for atmospheric Hg chemistry modeling include Br concentrations, stability and reactions of HgI, and speciation and photoreduction of HgII in aerosols and clouds.

Theoretical Study of the Monohydration of Mercury Compounds of Atmospheric Interest.

The structures, vibrational frequencies, and model IR spectra of the monohydrates of oxygenated mercury compounds (BrHgO, BrHgOH, BrHgOOH, BrHgNO2, BrHgONO, and HgOH) have been theoretically studied using the ωB97X-D/aug-cc-pVTZ level of theory. The ground state potential energy surface exhibits several stable structures of these monohydrates. The thermodynamic properties of the hydration reactions have been calculated at different levels of theory including DFT and coupled-cluster calculations DK-CCSD(T) with the ANO-RCC-Large basis sets. Standard enthalpies and Gibbs free energies of hydration were computed. The temperature dependence of ΔrG°(T) was evaluated for the most stable complexes over the temperature range 200-400 K. Thermodynamic data revealed that the highest fraction hydrated at 298 K and 100% relative humidity will be BrHgNO2-H2O at ∼5%.

On the work function of the surface Mo(0 0 1) and its temperature dependence: an ab initio molecular dynamics study.

Ab initio molecular dynamics simulations in NVT ensemble have been performed to investigate the finite temperature structure of the Mo(0 0 1) surface and its effect on work function (ϕ). In accord with previous experimental and theoretical work, our simulations predict that a termination with a stable reconstruction pattern is formed at T = 123 K. This pattern vanishes when temperature is increased to 423 K or 623 K and a disordered surface phase is formed whose time average corresponds to a bulk-like termination. Our results demonstrate that the surface relaxation is an important factor contributing to thermal variation of ϕ. At the lowest temperature, at which a stable reconstruction pattern is formed, the work function is found to increase by ∼0.23 eV compared to relaxed unreconstructed surface. The disappearance of stable reconstruction pattern at elevated temperatures leads to a decrease of ϕ by ∼0.07 eV. In contrast, the values computed for a non-reconstructing surface Mo(1 1 0) at T = 123 K, 423 K and 623 K are found to be nearly identical to the zero temperature value, which is a consequence of restricted atomic motion due to high packing density in this surface.

Theoretical Study of the Reactions of H Atoms with CH3I and CH2I2.

High level ab initio methods have been used to provide reliable kinetic data for the H + CH3I and H + CH2I2 gas-phase reactions. The (H, I)-abstraction and I-substitution reaction pathways were identified. The structures were determined on the potential energy surface at the MP2/aug-cc-pVTZ level of theory. The energetics was then refined using the coupled cluster theory. For the iodinated species, the spin-orbit coupling was calculated using the MRCI approach. The core valence and the scalar relativistic corrections were considered. Thermal rate constants were reported using the canonical transition-state theory (TST) and compared to computed values with the canonical variational transition-state theory (CVT) using the zero curvature tunneling (ZCT) and the small curvature tunneling (SCT) corrections over a wide temperature range (250-2500 K) to show the importance of quantum tunneling effects at low temperatures. They are given by the following expressions for the overall reactions using the CVT/SCT method: kH+CH3I( T) = 1.07 × 10-17 × T2.13 exp(2.68 (kJ mol-1)/ RT) and kH+CH2I2( T) = 5.73 × 10-21 × T2.97 exp(3.15 (kJ mol-1)/ RT). The I-abstraction is predicted to be the major pathway for both H + CH3I and H + CH2I2 reactions. The obtained kinetic parameters for the H + CH3I reaction are in excellent agreement with their experimental counterparts over the temperature range 300-750 K. On the basis of our calculated reaction enthalpies, a new evaluation of the standard enthalpy of formation at 298 K of CH2I and CHI2 has been provided. Obtained values are Δf H°298K (CH2I) = 219.5 kJ mol-1 and Δf H°298K(CHI2) = 296.3 kJ mol-1.

Reactivity of Hydrogen Peroxide with Br and I Atoms.

The reaction mechanisms of Br and I atoms with H2O2 have been investigated using DFT and high-level ab initio calculations. The H-abstraction and OH-abstraction channels were highlighted. The geometries of the stationary points were optimized at the B3LYP/aug-cc-pVTZ level of theory, and the energetics were recalculated with the coupled cluster theory. Spin-orbit coupling for each halogenated species was also explicitly computed by employing the MRCI level of theory. Thermochemistry for HOBr and HOI has been revised and updated standard enthalpies of formation at 298 K for HOBr and HOI are the following: ΔfH°298K(HOBr) = (-66.2 ± 4.6) kJ mol-1 and ΔfH°298K(HOI) = (-66.8 ± 4.7) kJ mol-1. The rate constants have been estimated using transition state theory (TST), canonical variational transition state theory (CVT), and CVT with small curvature tunneling (CVT/SCT) over a wide temperature range (250-2500 K). For the direct abstraction mechanism, the overall rate constant at 300 K was predicted to be 2.58 × 10-16 and 7.42 × 10-25 cm3 molecule-1s-1 for the Br + H2O2 and I + H2O2 reactions, respectively. The modified Arrhenius parameters have been estimated for the overall reactions: kBr+H2O2(T) = 4.80 × 10-26 T4.31 exp(-5.51 (kJ mol-1)/RT) and kI+H2O2(T) = 3.41 × 10-23 T3.29 exp(-56.32 (kJ mol-1)/RT).

MS-CASPT2 study of the ground and low lying states of CsH.

Correlated ab initio methods [CASPT2 and R-CCSD(T)] in conjunction with the ANO-RCC basis sets in large contraction were used to calculate potential energy curves (PECs) of the ground and excited electronic states of CsH+ (doublets and quartets) with the inclusion of the scalar relativistic effects and spin-orbit interaction. The ground X2Σ+ state is a rather fragile van der Waals molecular ion. The binding energy of this X2Σ+ state provided by both computational methods is estimated to be 0.02-0.04 eV, and is compared with the reported experimental binding energy (0.51-0.77 eV). This large binding energy can be attributed to the A2Σ+ state, and can thus explain the apparent disagreement between theory and experiment. The spectroscopic constants of all bound states were calculated from the PECs and compared with previous published data for X2Σ+ and A2Σ+ states. Graphical abstract Low-lying Ω states of cesium hydride cation.

Caesium hydride: MS-CASPT2 potential energy curves and A1Σ+→X1Σ+ absorption/emission spectroscopy.

Correlated ab initio methods (CASPT2 and CCSD(T)) in conjunction with the ANO-RCC basis sets were used to calculate potential energy curves (PECs) of the ground, valence, and Rydberg electronic states of CsH with the inclusion of the scalar relativistic effects. The spectroscopic constants of bound states were calculated from the PECs and compared with previous theoretical and/or available experimental data. Absorption and emission spectra arising from the transition between X1Σ+ and A1Σ+ states were modelled using vibrational and rotational energy levels and corresponding nuclear wave functions obtained via the direct numerical integration of one-dimensional rovibrational Schrödinger equation in the CASPT2/ANO-RCC electronic potentials. The anharmonic shape of the A1Σ+ potential and the shape of the pertinent vibrational wave functions have an interesting impact on the final shape of the spectrum and result in the complicated fine structure of individual emission bands.

A Density Functional Theory and ab Initio Investigation of the Oxidation Reaction of CO by IO Radicals.

To get an insight into the possible reactivity between iodine oxides and CO, a first step was to study the thermochemical properties and kinetic parameters of the reaction between IO and CO using theoretical chemistry tools. All stationary points involved were optimized using the Becke's three-parameter hybrid exchange functional coupled with the Lee-Yang-Parr nonlocal correlation functional (B3LYP) and the Møller-Plesset second-order perturbation theory (MP2). Single-point energy calculations were performed using the coupled cluster theory with the iterative inclusion of singles and doubles and the perturbative estimation for triple excitations (CCSD(T)) and the aug-cc-pVnZ (n = T, Q, and 5) basis sets on geometries previously optimized at the aug-cc-pVTZ level. The energetics was then recalculated using the one-component DK-CCSD(T) approach with the relativistic ANO basis sets. The spin-orbit coupling for the iodine containing species was calculated a posteriori using the restricted active space state interaction method in conjunction with the multiconfigurational perturbation theory (CASPT2/RASSI) employing the complete active space (CASSCF) wave function as the reference. The CCSD(T) energies were also corrected for BSSE for molecular complexes and refined with the extrapolation to CBS limit while the DK-CCSD(T) values were refined with the extrapolation to FCI. The exploration of the potential energy surface revealed a two-steps mechanism with a trans and a cis pathway. The rate constants for the direct and complex mechanism were computed as a function of temperature (250-2500 K) using the canonical transition state theory. The three-parameter Arrhenius expressions obtained for the direct and indirect mechanism at the DK-CCSD(T)-cf level of theory is 1.49 × 10(-17) × T(1.77) exp(-47.4 (kJ mol(-1))/RT).

Reactivity of CHI3 with OH radicals: X-abstraction reaction pathways (X = H, I), atmospheric chemistry, and nuclear safety.

The X-abstraction (X = H, I) pathways in the reaction of CHI3 with OH radical, a possible iodoform removal process relevant to the Earth's atmosphere and conditions prevailing in the case of a nuclear accident, have been studied applying highly correlated ab initio quantum chemistry methods and canonical transition-state theory to obtain reaction energy profiles and rate constants. Geometry optimizations of reactants, products, molecular complexes, and transition states determined at the MP2/cc-pVTZ level of theory have been followed by DK-CCSD(T)/ANO-RCC single-point energy calculations. Further improvement of electronic energies has been achieved by applying spin-orbit coupling, corrections toward full configuration interaction, vibration contributions, and tunneling corrections. Calculated reaction enthalpies at 0 K are -108.2 and -5.1 kJ mol(-1) for the H- and I-abstraction pathways, respectively; the strongly exothermic H-abstraction pathway is energetically favored over the modestly exothermic I-abstraction one. The overall rate constant at 298 K based on our ab initio calculations is 4.90 × 10(-11) cm(3) molecule(-1) s(-1), with the I-abstraction pathway being the major channel over the temperature range of 250-2000 K. The CHI3 atmospheric lifetime with respect to the removal reaction with OH radical is predicted to be about 6 h, very short compared to that of other halomethanes.

Microhydration of caesium compounds: Cs, CsOH, CsI and Cs2I2 complexes with one to three H2O molecules of nuclear safety interest.

Structure and thermodynamic properties (standard enthalpies of formation and Gibbs free energies) of hydrated caesium species of nuclear safety interest, Cs, CsOH, CsI and its dimer Cs2I2, with one up to three water molecules, are calculated to assess their possible existence in severe accident occurring to a pressurized water reactor. The calculations were performed using the coupled cluster theory including single, double and non-iterative triple substitutions (CCSD(T)) in conjunction with the basis sets (ANO-RCC) developed for scalar relativistic calculations. The second-order spin-free Douglas-Kroll-Hess Hamiltonian was used to account for the scalar relativistic effects. Thermodynamic properties obtained by these correlated ab initio calculations (entropies and thermal capacities at constant pressure as a function of temperature) are used in nuclear accident simulations using ASTEC/SOPHAEROS software. Interaction energies, standard enthalpies and Gibbs free energies of successive water molecules addition determine the ordering of the complexes. CsOH forms the most hydrated stable complexes followed by CsI, Cs2I2, and Cs. CsOH still exists in steam atmosphere even at quite high temperature, up to around 1100 K.

On the way to the highest coordination number in the planar metal-centred aromatic Ta©B10- cluster: evolution of the structures of TaB(n)- (n = 3-8).

The structures and chemical bonding of TaB(n)(-) (n = 3-8) clusters are investigated systematically to elucidate the formation of the planar metal-centred aromatic borometallic cluster, Ta©B10(-) (the © sign is used to designate the central position of the doped atom in monocyclic structures in M©B(n)-type planar clusters), which was found previously to have the highest coordination number for a metal atom in a planar geometry. Photoelectron spectroscopy is combined with ab initio calculations to determine the global minima of the TaB(n)(-) clusters. We find that from TaB3(-) to TaB5(-) the boron atoms nucleate around the central Ta atom to form fan-like structures. A structural transition occurs at TaB6(-), which is found to have a hexagonal structure, but with a boron atom in the centre and the Ta atom on the periphery. TaB7(-) is shown to have a three-dimensional boat-like structure, which can be viewed as a Ta atom coordinated to an elongated B7 cluster from above. The global minimum of the TaB8(-) cluster is found to be pyramidal with the Ta atom interacting with a B8 monocyclic ring. Starting from this structure, additional boron atoms simply enlarge the boron ring to form the slightly pyramidal TaB9(-) cluster and eventually the perfectly planar Ta-centred B10-ring aromatic cluster, Ta©B10(-). It is shown that boron atoms do not nucleate smoothly around a Ta atom on the way to the decacoordinated Ta©B10 (-) molecular wheel, but rather the competition between B-B interactions and Ta-B interactions determines the most stable structures of the smaller TaB(n)(-) (n = 3-8) clusters.

Structural changes in the series of boron-carbon mixed clusters C(x)B(10-x)- (x = 3-10) upon substitution of boron by carbon.

We report a theoretical investigation on the ten-atom boron-carbon mixed clusters C(x)B(10-x)(-) (x = 3-10), revealing a molecular wheel to monocyclic ring and linear species structural change as a function of x upon increasing the number of carbon atoms in the studied series. The unbiased searches for the global minimum structures of the clusters with x ranging from 3 to 9 were conducted using the Coalescence Kick program for different spin multiplicities. Subsequent geometry optimizations with follow-up frequency calculations at the hybrid density functional B3LYP∕6-311+G(d) level of theory along with the single point coupled-cluster calculations (UCCSD(T)/aug-cc-pVTZ//B3LYP/6-311+G(d) and RCCSD(T)/aug-cc-pVTZ//B3LYP/6-311+G(d)) revealed that the C3B7(-) and C4B6(-) clusters possess planar distorted wheel-type structures with a single inner boron atom, similar to the recently reported CB9(-) and C2B8(-). Going from C5B5(-) to C9B(-) inclusive, monocyclic and ring-like structures are observed as the most stable ones on the PES. The first linear species in the presented series is found for the C10(-) cluster, which is almost isoenergetic with the one possessing a monocyclic geometry. The classical 2c-2e σ bonds are responsible for the peripheral bonding in both carbon- and boron-rich clusters, whereas multicenter σ bonding (nc-2e bonds with n > 2) on the inner fragments in boron-rich clusters is found to be the effective tool to describe their chemical bonding nature. It was shown that the structural transitions in the C(x)B(10-x)(-) series occur in part due to the preference of carbon to form localized bonds, which are found on the periphery of the clusters. Chemical bonding picture of C10(-) is explained on the basis of the geometrical structures of the C10 and C10(2-) clusters and their chemical bonding analyses.

Photoelectron spectroscopy and ab initio study of boron-carbon mixed clusters: CB9- and C2B8-.

We performed a joint photoelectron spectroscopy and ab initio study of two carbon-doped boron clusters, CB(9)(-) and C(2)B(8)(-). Unbiased computational searches revealed similar global minimum structures for both clusters. The comparison of the experimentally observed and theoretically calculated vertical detachment energies revealed that only the global minimum structure is responsible for the experimental spectra of CB(9)(-), whereas the two lowest-lying isomers of C(2)B(8)(-) contribute to the experimental spectra. The planar "distorted wheel" type structures with a single inner boron atom found for CB(9)(-) and C(2)B(8)(-) are different from the quasi-planar structure of B(10)(-), which consists of two inner atoms and eight peripheral boron atoms. The adaptive natural density partitioning chemical bonding analysis revealed that CB(9)(-) and C(2)B(8) clusters exhibit π aromaticity and σ antiaromaticity, which is consistent with their planar distorted structures.

Ab initio calculations of molecular properties of low-lying electronic states of 2-cyclopenten-1-one--link with biological activity.

Geometries, vibrational frequencies, vertical and adiabatic excitation energies, dipole moments and dipole polarizabilities of the ground and the three lowest electronic excited states, S1(n, π*), T1(n, π*), and T2(π, π*) of the 2-cyclopenten-1-one molecule (2CP) were calculated at the CCSD and CCSD(T) levels of approximation. Our results indicate that two triplets T1(n, π*) and T2(π, π*) are lying very close each to other, while the singlet S1(n, π*) is well above them. There are dramatic changes in dipole moments for (n, π*) excited states in respect to the ground state. On the other hand the T2(π, π*) state has a similar dipole moment as the ground state. These changes can be interpreted within the MO picture using electrostatic potential maps and changes in model IR spectra. Our CCSD(T) dipole moment data for the ground state and almost isoenergetic triplets T1(n, π*) and T2(π, π*) are 1.469 a.u., 0.551 a.u., and 1.124 a.u., respectively. Dipole polarizabilities of investigated excited states are much less affected by electron excitations than dipole moments. These are the first dipole moment and polarizability data of 2CP in the literature. The changes of molecular properties upon excitation to S1(n, π*) and T1(n, π*) correlate with the experimental data on the biological activity of 2CP related to the α, ß-unsaturated carbonyl group.

Toward more efficient CCSD(T) calculations of intermolecular interactions in model hydrogen-bonded and stacked dimers.

Interaction energies of the model H-bonded complexes, the formamide and formamidine dimers, as well as the stacked formaldehyde and ethylene dimers are calculated by the coupled cluster CCSD(T) method. These systems serve as a model for H-bonded and stacking interactions, typical in molecules participating in biological systems. We use the optimized virtual orbital space (OVOS) technique, by which the dimension of the space of virtual orbitals in coupled cluster CCSD(T) calculations can be significantly reduced. We demonstrate that when the space of virtual orbitals is reduced to 50% of the full space, which means reducing computational demands by 1 order of magnitude, the interaction energies for both H-bonded and stacked dimers are affected by no more than 0.1 kcal/mol. This error is much smaller than the error when interaction energies are calculated using limited basis sets.